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ABSTRACT 


An  investigation  into  the  detection  capability  of  a prototype 
seismic  network  composed  of  a 32  station  subset  of  the  approximately  150 
stations  used  by  the  ISM  experiment  has  been  completed.  This  investiga- 
tion was  threefold  in  nature.  The  first  investigation  directly  compared 
the  lists  of  events  produced  by  the  32  station  network  experiment,  the 
USGS-PDE,  and  the  ISM  experiment  for  the  ISM  time  period,  20  February 
1972  to  19  March  1972.  The  number  of  stations  used,  the  quality  of 
available  data  and  events  acceptance  criteria  were  not  the  same  for  all 
three  lists.  The  comparisons  determined  the  effects  the  different  input 
data  and  acceptance  criteria  had  on  the  number  of  events  reported  and 
the  relative  location  accuracy  of  each.  Secondly,  an  investigation  to 
directly  determine  the  effect  a careful  reading  of  seismograms  had  on 
the  performance  of  individual  stations  and/or  small  networks  was  conducted. 
To  evaluate  the  improvement  in  performance  for  single  stations,  2 stations 
were  selected  and  reread.  The  reread  detect'  s were  compared  to  the 
detections  that  were  reported  to  the  USGS  for  all  ISM  events  that  were 
within  the  distance  range  between  20°  and  90°  from  the  station.  These 
comparisons  help  verify  that  a significant  number  valid  and  clear 
arrivals  are  not  reported  to  USGS.  To  evaluate  the  impact  on  network 
performance  21  stations  of  the  32  station  network  were  reread  for  a 
small  number  of  events  to  determine  the  extent  to  which  the  rereading  of 
the  seismograms  could  change  the  detection  and  locating  capabilities  of 
a small  network.  The  final  investigation,  to  compliment  the  experimental 
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evaluations  was  to  program  a simple  statistical  model  to  predict  the 
probability  of  an  N station  detection  from  an  M station  network  with  a 
known  average  station  probability.  The  data  from  the  32  station  set 
were  used  to  estimate  the  single  station  detection  probability  for 
events  with  m^  magnitudes  of  4.6  but  less  than  4.7  within  a distance  of 
90°  or  less  from  the  station.  This  probability  was  entered  into  the 
model  and  a network  detection  probability  table  was  generated. 
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Introduction 


In  January  1972  a group  of  seismologists  from  several  countries  met 
in  Cambridge,  Massachusetts,  to  discuss  problems  related  to  the  seismic 
detection  and  identification  of  underground  nuclear  tests.  It  was 
decided  that  a cooperative  attempt  be  made  to  better  assess  the  capabil- 
ities of  deployed  seismic  instrumentation  for  monitoring  world  wide 
seismic  activity.  The  period  from  20  February  to  19  March  1972  (the 
International  Seismic  Month-ISM)^ was  selected  for  this  experiment 
with  the  Lincoln  Laboratory  acting  as  a data  center. 

A list  of  996  seismic  events  were  reported  by  Lincoln  Laboratory^ 
as  a result  of  the  experiment.  A secondary  result  of  the  ISM  experiment 
has  been  the  evaluation  of  the  capability  of  a small  network  of  stations 


to  monitor  global  seismic  activity.  Thirty-two  seismic  stations  from 
the  approximately  150  used  in  the  ISM  experiment  were  selected  to  be 
evaluated  as  a small  network.  The  32  stations  were  chosen  ad  hoc  to 
give  good  global  coverage  while  using  stations  which  contributed  data  to 
a large  number  of  ISM  events.  The  results  of  an  in  depth  comparison  of 
the  list  of  events  produced  using  the  network,  the  ISM  experiment  list, 
and  the  USGS  PDE  list  (United  States  Geological  Survey  - Preliminary 
Determination  of  Epicenters)  are  given  below.  To  better  evaluate  the 
detection  capability  of  a small  network  of  32  stations,  two  additional 
interrelated  experiments  were  also  completed  and  are  reported  below. 

One  of  these  experiments  relates  to  the  improvement  which  can  be  achieved 
in  the  event  detection  capability  of  individual  stations  which  report  to 
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the  USGS  by  a more  exhaustive  reading  of  the  seismograms.  The  other  was 
a simple  theoretical  study  of  a small  network’s  probability  of  detecting 
and  locating  an  event. 

The  three  lists  discussed  in  this  report  are  composed  of  events 
occurring  during  the  ISM  experiment  time  period,  20  February  1972  through 
19  March  1972.  The  methods  used  and  the  parameters  needed  for  an  event 
acceptance  for  each  of  these  lists  are  as  follows: 


USGS  PDE  List 

The  USGS  uses  arrival  times  from  approximately  600  seismograph 

stations  to  prepare  its  event  list.  The  criteria  required  for 

event  acceptance  are  that  at  least  five  station  initial  arrival 

times  (P  or  PKP)  are  associated  from  at  least  two  different  azimuths. 

For  near  regional  events,  phases  P , P , S , and  S can  be  accepted. 

n g n g 

Not  all  events  meeting  these  acceptance  requirements  are  necessarily 
reported.  An  event  can  be  rejected  at  the  discretion  of  a USGS 
analyst  who  reviews  the  data. 


ISM  List 

The  ISM  experiment  used  data  reported  to  USGS  only  for  stations 
which  reported  at  least  50  times  during  the  ISM  period.  This 
requirement  eliminated  a large  majority  of  the  stations  that  report 
to  the  USGS,  leaving  approximately  150  stations  used  for  the  ISM 
experiment.  All  of  the  stations  used  in  the  computation  of  events 
by  the  ISM  experiment,  with  the  exception  of  stations  CHG  and  YKA, 
were  contributers  to  the  PDE  computations,  but  not  necessarily  at 
the  same  detection  level.  The  detections  which  the  USGS  used  from 
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the  two  large  arrays,  NAO  and  LAO,  were  from  the  published  processed 
event  bulletins.  The  ISM  experiment  used  these  plus  preliminary 
detection  logs  which  contained  many  arrivals  which  do  not  appear  in 
the  site  bulletins.  Data  from  the  array  YKA  were  not  reported  to 
the  USGS  but  all  detections  were  processed  for  use  in  the  ISM 
experiment.  Swedish  seismologists  reprocessed  HFS  data  and  reported 
many  more  arrivals  to  the  ISM  experiment  than  were  reported  to  the 
USGS.  All  Canadian  stations  were  reread  by  personnel  of  the  Seismo- 
logical  section  of  the  Canadian  Department  of  Energy,  Mines  and 
Resources,  and  these  were  used  in  the  ISM  experiment.  The  number 
of  reread  detections  from  the  Canadian  stations  greatly  exceeded 
the  number  of  detections  sent  to  the  USGS.  The  stations  KBL  and 
CHG  were  reread  from  film  chips  by  an  analyst  within  the  Lincoln 
Laboratory  group.  The  detections  for  the  rest  of  the  stations  used 
in  the  ISM  experiment  were  exactly  the  same  as  were  reported  to  the 
USGS. 

The  ISM  experiment  event  acceptance  criteria  differ  from  the 
PDE  requirements  in  that  only  4 detections,  instead  of  are 
required  and  they  need  not  be  just  initial  arrivals  of  phases  P or 
PKP.  The  ISM  locate  program  accepts  phases  P,  PP,  PcP,  ScP,  pP,  3 
branches  of  PKP,  3 branches  of  PKKP,  SKP  and  PfPf  detections.  Also 
listed  as  acceptable  detections  are  array  computed  dt/dA  and  array 
computed  azimuths.  The  ISM  experiment  also  requires  that  at  least 
one  of  the  associated  stations  have  a distance ^ 3°  from  the  epicenter. 
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Fig.  1.  Geographic  distribution  of  prototype  network 
stations.  Not  shown  on  this  map  is  station  SPA  (South 
Pole,  Antartica). 


The  final  acceptance  of  events  listed  on  the  ISM  list  was  at  the 
discretion  of  Lincoln  Laboratory  analysts  after  reviewing  the  event 
location  parameters. 

32  Station  List 

The  seismic  stations  selected  to  make  up  the  32  station  network 
are  listed  on  Table  1 and  a map  showing  their  geographic  distribution 
is  shown  on  Figure  1.  This  32  station  network  is  a subset  of  the 
ISM  network  and  the  event  acceptance  criteria  were  the  same  as  for 
the  ISM  list.  For  a given  station  the  same  detections  were  available 
as  for  the  ISM  experiment.  No  new  detections  were  added. 
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Table  1 

32  Station  Network 


LAO 

YKA 

UBO 

NAO 

HFS 

MBC 

KBL 

ASP 

MAT 

COL 

CHG 

PNS 

CTA 

BLC 

NUB 

TUC 

FBC 

GBA 

SPA 

PMG 

KIC 

CLL 

SSF 

BDF 

BNG 

QUE 

SHI 

BAG 

BUL 

AFI 

EZN 

SHL 


LASA  Array,  Montana 

Yellowknife  Array,  N.  W.  Terr.  Canada 
Unita  Basin  Observatory,  Utah 
NORSAR  Array,  Norway 
Hagfors  Array,  Sweden 
Mould  Bay,  N.  W.  Terr.  Canada 
Kabul,  Afghanistan  (WWSSN) 

Alice  Springs,  N.  Terr.  Australia 
Matsushiro,  Honshu,  Japan  (WWSSN) 

College  Outpost,  Alaska  (WWSSN) 

Chengmai,  Thialand  (WWSSN) 

Penas,  Bolivia 

Charters  Towers,  Queensland,  Australia  (WWSSN) 
Baker  Lake,  N.  W.  Terr.  Canada 
Nurmi jarvi,  Finland  (WWSSN) 

Tucson,  Arizona  (WWSSN) 

Frobisher  Bay,  Canada 
Gauribidanur  Array,  India 
South  Pole,  Antarctica  (WWSSN) 

Port  Moresley,  New  Guinea  (WWSSN) 

Koson  Boka,  Ivory  Coast 
Collmberg,  E.  Germany 
Saint  Saulge,  France 
Brasilia  Array,  Brazil  (WWSSN) 

Bangui,  Central  Africa  Rep. 

Quetta,  Pakistan  (WWSSN) 

Shiraz,  Iran  (WWSSN) 

Baguio  City,  Philippines  (WWSSN) 

Bulawayo,  Rhodesia  (WWSSN) 

Afiamalu,  Samoa  (WWSSN) 

Ezine,  Turkey 
Shillong,  India  (WWSSN) 
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Event  List  Comparisons 


The  PDE  list  contained  353  events  for  the  ISM  time  period.  Only 

two  of  these  did  not  appear  in  the  ISM  list.  The  two  events  not  reported 

in  the  ISM  list  were  a small  Alaskan  event  with  no  reporting  station 

greater  than  3°  from  the  epicenter  and  a small  South  African  event  which 

required  the  use  of  phase  P , an  option  which  the  ISM  locate  program  did 

6 

not  have. 

If  the  USGS  had  used  only  the  data  submitted  to  them  from  stations 
that  were  used  by  the  ISM  experiment,  then  3^6  events  would  have  met  the 
PDE  event  requirements.  The  seven  events  that  would  not  have  been 
reported,  were  small  events  detected  and  located  by  a network  of  close- 
in  stations  to  the  epicenter,  in  most  cases  <10°.  Four  of  these  seven 
events  had  no  magntiude  computed  while  the  other  three  had  magnitudes  of 
3.8,  3.8,  and  U.9,  with  the  U.9  m^  being  a single  station  value  for  a 
station  10°  from  the  epicenter. 

We  don't  know  how  many  events  the  USGS  would  have  obtained  if  they 
had  limited  themselves  to  the  stations  used  by  the  ISM  but  had  used  all 
phase  and  velocity  information  reported  to  them.  However,  a check  of 
three  months  of  data  reported  and  entered  into  the  PDE  hypocenter 

7 

system  indicates  that  a monthly  average  of  11,303  secondary  phases, 

17 1+1  array  computed  velocities  and  17^+1  array  computed  azimuths  could 
have  been  added  to  the  initial  arrival  detection  base  used  by  USGS  for 
event  association  and  acceptance. 
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Of  the  353  PDE  events,  308  could  be  located  using  data  from  the  32 
station  subset  and  using  ISM  event  acceptance  criteria.  A histogram  of 
the  ISM  AVEMB  (average  magnitude)  of  the  ^3  events  reported  by  the 
PDE  list  but  not  reported  by  the  32  station  list  is  shown  in  Figure  2. 

The  shaded  portions  of  Figure  2 are  events  which  had  one  or  less  individual 
station  m^ 1 s averaged  into  AVEMB.  One  can  note  that  only  2 events  with 
an  greater  than  b . 5 were  not  reported  by  the  32  station  list,  and 
each  of  these  events  have  a questionable  single  station  AVEMB  value. 

The  ISM  list  is  composed  of  996  events  which  occurred  during  the 
ISM  period.  Of  these  996  events,  787  events  satisfy  the  USGS  event 
acceptance  criteria.  This  increase  of  ^3^-  events  over  the  353  reported 
on  the  PDE  list  is  due  to  the  improved  reporting  of  the  « 25  reread 
stations.  A histogram  of  the  ISM  AVEMB  distribution  of  the  6^3  events 
not  reported  by  the  USGS  is  shown  in  Figure  3.  The  shaded  portion 
denotes  the  events  that  had  data  from  less  than  2 stations  used  to 
calculate  the  AVEMB.  The  majority  of  the  events  greater  than  U.8  m^  are 
events  with  single  station  values  which  are  always  somewhat  question- 
able . 

The  32  station  network,  using  ISM  event  acceptance  criteria  reported 
8l8  of  the  996  events  reported  on  the  ISM  list.  Of  these  8l8  events  U60 
would  have  been  reported  if  the  event  acceptance  requirements  of  the 
USGS  had  been  used.  This  increase  of  107  events  over  the  353  events 
reported  by  the  entire  USGS  network  on  the  PDE  list  is  due  to  improved 
detection  from  only  b arrays  and  5 single  stations.  The  USGS,  using  the 
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Fig.  2.  Magnitude  distribution  of  PDE  events  not 
reported  by  the  prototype  network. 
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Fig.  3.  Magnitude  distribution  of  ISM  reported  events 
not  reported  by  PDF. 
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detections  from  these  32  stations  which  were  reported  to  them  would  have 
reported  only  25^  events.  A histogram  of  the  ISM  AVEMB  distribution  of 
the  508  32  station  list  events  not  reported  by  the  USGS  is  shown  as 
Figure  h with  the  shaded  portion  being  events  with  the  less  reliable 
single  station  average  m^  values.  The  number  of  events  reported  on  the 
32  station  list  was  177  events  less  than  the  number  of  events  on  the  ISM 
list.  A histogram  of  the  ISM  AVEMB  distribution  of  these  events  is 
shown  as  Figure  5.  Only  23  of  the  177  events  have  an  average  m^  that 
used  2 or  more  station  ni^’s  in  its  average.  Many  of  the  177  events  of 
the  ISM  list  and  the  U3  events  of  the  PDE  list,  not  reported  by  the  32 
station  list,  were  small  events  requiring  a more  compact  geographically 
distributed  network  than  the  32  station  network  to  meet  event  acceptance 
criteria  of  this  experiment. 

To  determine  the  accuracy  of  the  epicenter  locations  reported  on 
the  ISM  list  of  events  and  the  32  station  list  of  events,  the  following 
comparisons  were  made.  A direct  comparison  of  the  epicenter  location 
parameters  of  the  351  events  common  to  both  the  PDE  list  and  the  ISM 
list  was  made  which  indicated  that  99%  of  the  ISM  event  epicenters  were 
within  1°  of  the  PDE  epicenters.  As  was  previously  stated,  787  of  the 
996  ISM  events  would  have  satisfied  the  USGS  PDE  event  acceptance  criteria. 
Assuming  that  the  extra  U36  events  would  therefore  have  locations  as 
good  as  those  of  the  351  directly  compared  events  we  conclude  that  these 
787  ISM  events  would  have  epicenter  locations  within  1°  of  the  PDE 
epicenters.  Using  this  expanded  base  of  reliable  epicenters,  a direct 
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Fig.  U.  Magnitude  distribution  of  prototype  network 
events  not  reported  by  PDE. 
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Fig.  5*  Magnitude  distribution  of  ISM  events  not 
reported  by  the  prototype  network. 
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comparison  of  the  U60  32-station  event  epicenters  to  common  events  of 
the  787  ISM  events  was  made.  Figure  6 shows  the  epicenter  errors  obtained 
from  this  comparison,  which  indicate  that  96%  of  these  events  have 
locations  within  1°  of  the  ISM  epicenters. 

The  following  conclusion  can  be  made  from  these  event  list  comparison 
statistics.  (l)  A large  network  of  stations  approximately  1/k  the  size 
of  the  USGS  network,  such  as  the  ISM  network,  using  improved  detection 
methods  from  less  than  20%  of  its  stations  and  using  different  event 
acceptance  criteria  can  vastly  increase  the  number  of  reported  events  in 
the  magnitude  range  4.5  ni^  (Figure  3)  without  appreciable  deterioration 
in  epicenter  location  accuracy.  (2)  A small  network  of  stations,  such 
as  the  32  station  network,  with  approximately  30%  of  these  stations 
reporting  reread  detections  and  using  the  event  acceptance  criteria 
required  by  the  USGS,  can  not  only  increase  the  number  of  reported 
events  but  also  have  epicenters  that  are  accurate  to  within  approximately 
1°  of  the  PDE  epicenters. 
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Performance  Improvements  Resulting  from  Careful  Rereading  of  Seismograms 


Two  stations  of  the  32  station  network  were  selected  to  demonstrate 
by  example  how  typical  stations  are  not  currently  used  to  full  capability. 
Stations  MBC  and  KBL  were  chosen  as  stations  whose  data  was  reported  to 
both  the  USGS  and  the  ISM  experiment.  The  detections  from  MBC  reported 
to  the  ISM  were  reread  by  members  of  the  Canadian  Department  of  Energy, 
Mines  and  Resources,  while  those  from  KBL  were  reread  by  an  analyst  of 
the  Lincoln  Laboratory  group.  A comparison  was  made  of  the  detections 
reported  to  the  USGS  and  to  the  ISM  experiment  for  ISM  events  20°  ^ A < 

90°  from  each  of  these  stations.  We  also  restricted  the  study  to  events 
with  an  ISM  average  magnitude  determined  by  at  least  two  stations  to 
eliminate  confusion  resulting  from  poor  m^  values.  Figure  7 shows  the 
ISM  magnitude  distribution  of  all  ISM  events  within  this  distance  range 
for  station  MBC.  The  dotted  shading  of  this  histogram  denotes  the 
events  for  which  detections  were  sent  from  the  station  to  the  USGS  for 
their  use.  The  cross  line  shading  indicates  additional  detections 
obtained  by  the  more  thorough  reading  of  MBC  seismograms  which  was  done 
for  the  ISM.  For  magnitudes  of  <5*0  m^,  one  notes  a large  increase  in 
the  number  of  associated  detections  obtained  by  rereading.  Figure  8 
shows  the  associated  detections  of  Figure  7 converted  into  cumulative 
percentage  of  ISM  events.  The  dips  observed  on  the  cumulative  curves 
for  magnitude  greater  than  5.0  in^  are  due  to  a few  missed  events  and  are 
probably  an  artifact  resulting  from  the  small  sample  size  for  these 
larger  magnitudes.  This  artifact  makes  it  difficult  to  precisely 
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Fig.  7.  Magnitude  distribution  of  events  20°  A<90° 
from  MBC. 
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Fig.  8.  Cumulative  percentage  of  ISM  events  for  MBC 
detections . 
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compare  detection  levels.  Nevertheless  it  seems  clear  that  both  the  50% 
and  90%  detection  thresholds  are  at  least  1.0  m^  units  lower  for  the 
reread  data. 

The  same  comparison  was  made  for  station  KBL,  which  according  to 

g 

NEIS  (National  Earthquake  Information  Service)  of  the  USGS  was  one  of 
the  better  reporting  stations  for  PDE  computations.  The  event  to  station 
distance  restraints  and  the  number  of  station  m^  values  required  to 
compute  average  m^  were  the  same  for  KBL  as  was  stated  for  MBC.  Figure 
9 shows  this  comparison.  The  improvement  in  the  number  of  events  associ- 
ated with  the  reread  detections  is  significant  for  m^’s  < ^ . 5 • Figure  10 
shows  the  cumulative  detection  curves  for  KBL.  The  90%  cumulative 
detection  threshold  did  not  change  but  an  improvement  of  ^ 0.8  m^  units 
at  the  50%  level  is  clear.  Although  some  improvement  was  obtained  by 
rereading  the  seismograms  from  station  KBL,  it  is  not  as  great  as  that 
obtained  by  rereading  the  seismograms  from  station  MBC. 

Now  suppose  that  all  of  the  non-array  stations  in  the  32  station 
network  had  been  read  as  carefully  as  possible.  What  would  be  the 
capability  of  such  a network?  To  address  this  question  we  chose  limit 
ourselves  to  events  with  magnitudes  U.6  £ < U.7  (the  90%  cumulative 

detection  threshold  of  the  ISM  experiment^).  There  were  38  events  in 
this  magnitude  range  on  the  ISM  list.  Of  these  only  2h  were  used  on  the 
USGS  PDE  list.  Table  2 shows  individual  station  detection  statistics 
for  these  38  events.  The  detection  percentages  given  are  based  only  on 
the  subset  of  events  within  90°  of  the  particular  station.  The  number 
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Fig.  10.  Cumulative  percentage  of  ISM  events  for 
KBL  detections. 
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of  events  within  90°  of  each  station  is  presented  in  column  1.  Column  2 
indicates  the  number  of  these  events  which  had  detections  reported  to 
the  USGS.  Column  3 is  just  the  ratio  of  column  2 to  column  1.  This 
ratio  is  our  estimate  of  the  probability  that  the  station  will  report  an 
arrival  to  the  USGS  for  an  event  in  the  U.6  m^  magnitude  range  which  is 
less  than  90°  from  the  station.  The  average  probability  is  0.35. 

Column  h is  the  same  as  column  2 but  for  the  detections  used  in  the  ISM 
experiment.  Column  5 is  our  estimate  detection  probability  based  upon 
data  available  for  the  ISM.  The  average  probability  in  this  case  is 
increased  to  O.UU.  The  average  station  probability  improvement  for  the 
ISM  is  caused  by  the  lower  detection  threshold  of  only  9 of  the  32 
stations  (denoted  by  * in  column  U).  For  this  experiment,  12  additional 
stations  of  the  32  station  network  were  reread  from  film  chips  making  a 
total  of  21  of  the  32  stations  being  reread.  Film  chips  were  not  avail- 
able to  us  for  the  remaining  11  stations  denoted  with  a +.  Column  6 is 
similar  to  2 and  h but  for  this  most  complete  visual  analysis  of  the 
available  seismograms.  Column  7 is  the  station  detection  probability 
estimated  from  1 and  6.  The  average  probability  is  now  .563  and  if  only 
the  21  stations  that  were  actually  reread  are  averaged  the  average 
probability  is  .665. 

An  attempt  was  made  to  locate  the  38  ISM  events  with  U.6  ^ m^  < U.7 
using  these  various  sets  of  detections  at  the  32  station  network.  We 
used  both  the  USGS  PDE  event  acceptance  criteria  and  the  ISM  experiment 
event  acceptance  criteria.  Tables  3 and  h show  the  results. 
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Table  2 


Station  Detection  Statistics 


1 

No.  Events 

2 

3 

4 

5 

6 

7 

Station 

4.6  < M,  <4.7 
b 

A < 90° 

Reported 
To  USGS 

USGS 

Station 

Probability 

Reported 
To  ISM 

ISM 
Station 
Prababil  ity 

This 

Experiment 

Final 
Stotian 
Prababil  ity 

LAO 

16 

11 

0.688 

15* 

0.939 

15* 

0.939 

YKA 

24 

0 

0.000 

14* 

0.583 

14* 

0.583 

UBO 

17 

12 

0.706 

12 

0.706 

1 2 1 

0.706 

NAO 

21 

15 

0.714 

19* 

0.905 

19* 

0.905 

HFS 

22 

9 

0.409 

14* 

0.636 

14* 

0.636 

MBC 

20 

7 

0.350 

15* 

0.750 

15* 

0.750 

KBL 

21 

11 

0.524 

13* 

0.619 

13* 

0.619 

ASP 

27 

19 

0.704 

19 

0.704 

19t 

0.704 

MAT 

32 

10 

0.313 

10 

0.313 

23* 

0.719 

COL 

30 

15 

0.500 

15 

0.500 

23* 

0.767 

CHG 

27 

0 

0.000 

13* 

0.481 

20* 

0.741 

PNS 

8 

1 

0.125 

1 

0.125 

it 

0.125 

CTA 

29 

17 

0.586 

17 

0.586 

1 7 1 

0.586 

BLC 

18 

5 

0.279 

9* 

0.500 

9t 

0.500 

NUR 

20 

10 

0.500 

10 

0.500 

13* 

0.650 

TUC 

17 

4 

0.235 

4 

0.235 

10* 

0.588 

FBC 

17 

2 

0.119 

10* 

0.588 

10* 

0.588 

GBA 

20 

8 

0.400 

8 

0.400 

8t 

0.400 

SPA 

16 

9 

0.563 

9 

0.563 

15* 

0.938 

PMG 

28 

10 

0.357 

10 

0.357 

16* 

0.571 

KIC 

10 

5 

0.500 

5 

0.500 

5t 

0.500 

CLL 

21 

6 

0.286 

6 

0.286 

6t 

0.286 

SSF 

16 

4 

0.250 

4 

0.250 

4t 

0.250 

BDF 

8 

3 

0.375 

3 

0.375 

3t 

0.375 

BNG 

9 

5 

0.556 

5 

0.556 

5t 

0.556 

QUE 

21 

6 

0.273 

6 

0.273 

12 

0.571 

SHI 

17 

1 

0.059 

1 

0.059 

9* 

0.591 

BAG 

31 

6 

0.194 

6 

0.194 

14* 

0.452 

BUL 

9 

2 

0.222 

2 

0.222 

8* 

0.899 

AFI 

31 

5 

0.161 

5 

0. 161 

9* 

0.290 

EZN 

17 

1 

0.059 

1 

0.059 

it 

0.059 

SHL 

25 

2 

0.080 

0.347 

2 

0.080 

0.436 

17* 

0.680 

0.563 

0.665 

Statianst 

*Rereod  Stations 


tFilm  Chips  Not  Availoble  At  Lincoln  Laboratory  far  Rereading 
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Table  3 


32  Station  Set  Detection  Capability 


32  Station 
Detection 
Set 


Average 

Station 

Probability 


No.  of  38  Events 
Located  Using 
PDE  Acceptance 
Criteria 


% No.  of  38  Events  % 
Located  with  ISM 
Acceptance 
Criteria 


PDE 

.347 

21 

55.3 

30 

78.9 

ISM 

.436 

29 

76.3 

36 

94.7 

Reread 

ISM 

.563 

36 

94.7 

38 

100.0 

Table  h shows  the  location  capability  achieved  using  only  the 
21  stations  reread  for  this  experiment. 


Table  k 

21  Station  Set  Detection  Capability 


Detection 

Set 


Average 

Station 

Probability 


No.  of  38  Events 
Located  Using 
PDE  Acceptance 
Criteria 


No.  of  38  Events 
Located  with  ISM 
Acceptance 
Criteria 


21  STA 

.665 

32 

84.2 

35 

92.1 

2k 


These  data  verify  again  that  careful  reading  and  processing  of  20- 
BO  well  placed  stations  can  generate  a significantly  larger  event  list 
than  obtained  by  current  USGS  practice  which  makes  use  of  many  more 
stations.  Our  data  establishes  this  for  h.6  ^ m^  < U.7  where  the  USGS 
reporting  probability  is  low  (2U/38  = .63)-  For  higher  magnitudes  the 
difference  would  decrease  since  the  USGS  reporting  probability  would 
increase  toward  1.0.  Probably  the  difference  would  be  larger  for  smaller 
events  if  we  exclude  events  in  the  USGS  list  which  are  based  upon  only 
local  data.  Of  course* acceptance  criteria  also  significantly  alters 
event  list  detection  thresholds,  especially  with  low  detection  probabilities 
for  individual  stations.  These  acceptance  criteria  also  affect  location 
accuracy  but  we  have  not  addressed  this  question  here. 

Examples  of  additional  detections  that  were  obtained  by  rereading 
film  chips  are  shown  in  Figure  11  for  QUE.  Station  QUE  was  chosen  not 
as  an  example  of  poor  reading  and  reporting  but  as  an  average  reporting 
station.  Additional  detections  for  6 events  of  the  38  with  reliable  m^ 
values  were  obtained  by  rereading  this  station.  These  reread  detections 
are  the  result  of  hindsight  picking  but  at  least  h of  these  events  are 
quite  obvious.  Three  additional  events  not  shown  in  this  figure  could 
have  possibly  been  reported  by  a well  trained  analyst. 
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0439 


2/20  ORIGIN  TIME  04  30  42  AVERAGE  mR  4.63  QUEA48  580  THEORETICAL  ARRIVAL  TIME  AT  QUE  04  39  24 


W***Sm **^**#» *+y*\^**+4*m 

*^w***m^-^**w*f*t+m*^*mi**K  * m iu*  tMm******^*^  i iwl  i n 

2/23  ORIGIN  TIME  11  36  01  AVERAGE  mB4  63  QUE  A 75  32°  THEORETICAL  ARRIVAL  TIME  AT  QUE  11  47  42 


2316 


o 


rnmjw^d**  M/fi  NW*f  &mfp?**4**^ 


2/23  ORIGIN  TIME  23  13  42  AVERAGE  mB  4.68  QUE  A 12.77°  THEORETICAL  ARRIVAL  TIME  AT  QUE  23  16  43 


324»*N«  V1  _ 

3/02  ORIGIN  TIME  03:16:09  AVE  RAGE  mB  4.62  QUEA44.680  THEORETICAL  ARRIVAL  TIME  AT  QUE  03 : 24 : 1 1 
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» '^v  ***+*  v'-r**+* 


3/08  ORIGIN  TIME  15:23:00  AVERAGE  mB  4.69  QUE  A 85  41°  THEORETICAL  ARRIVAL  TIME  AT  QUE  15 : 35:  35 


0232 

3/10  ORIGIN  TIME  02:22:38  AVERAGE  mB  4.61  QUEA61.59C  THEORETICAL  ARRIVAL  TIME  AT  QUE  02 :32 : 52 
* AVERAGE  mB  OBTAINED  FROM  ONLY  1 mB  OBSERVATION 

Fig.  11.  Additional  detections  obtained  rereading 
station  QUE  film  chips. 
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Network  Detection  Probability 


A simple  statistical  model  to  predict  the  probability  of  an  N 
station  detection  from  an  M station  network  with  a known  average  station 
probability  has  been  programmed  to  complement  the  experimental  evaluation 
of  small  networks.  The  formula  used  for  this  model  is  as  follows: 


P 


M 

E 

k=N 


k,  *M-k 
P (1-p) 


where  P is  the  probability  of  N or  more  detections  of  an  event  by  M 
detectors  if  the  detectors  are  independent  and  each  has  a probability  p 
of  detecting  recognizing  any  particular  event. 

Suppose  p = .630.  This  average  station  probability  is  the  average 
of  the  best  reported  1 6 stations  of  the  32  station  network  for  events 
with  U.6^  < 4.7  which  are  within  90°  of  the  station  and  is  the 

approximate  average  one  would  achieve  if  all  of  the  32  station  network 
stations  were  reread.  Using  the  model  with  this  average  station  prob- 
ability, a network  detection  probability  table  has  been  generated  and  is 
shown  as  Table  5«  An  example  of  the  use  of  this  table  is  as  follows. 
Suppose  that  half  of  the  32  stations  are  within  90°  of  a particular 
event.  Then  the  probability  of  detection  by  5 stations  or  more  is 
0.998.  This  drops  to  0.907  for  8 or  more  and  0.625  for  10  or  more.  Of 
course  the  model  could  generate  such  a table  for  any  average  station 
probability.  The  particular  table  shown  indicates  that  a 32  station 
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network  should  have  a global  detection  capability  of  essentially  100$  at 
b.6,  which  is  the  estimated  90$  incremental  detection  threshold  for 
the  ISM  list. 
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Network  Detection  Probability 
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Conclusions  and  Discussions 


A network  consisting  of  a small  number  ( 2-3)  of  arrays  and  (20-30) 
traditional  seismic  stations  could  be  used  to  routinely  produce  a global 
seismic  bulletin  (excluding  local  events)  better  than  those  now  produced 
on  a routine  basis  by  USGS  or  the  large  arrays. 

A list  of  events  produced  using  data  from  this  small  network  would 
have  better  hypoc enter  control  than  do  lists  produced  by  the  individual 
seismic  arrays.  The  small  network  could  be  used  to  locate  considerably 
more  events  than  the  current  PDE  operation  does  without  sacrificing 
hypoc enter  accuracy. 

A most  significant  step  towards  achieving  the  potential  capability 
of  a small  network  would  be  simply  to  arrange  for  uniform  and  careful 
reading  of  seismograms  from  a few  tens  of  good  stations  reasonably 
distributed  about  the  globe.  Some  potential  advantages  of  the  better 
quality  control  on  the  data  submitted  from  stations  of  this  small  network, 
other  than  the  increased  number  of  detections,  would  be  more  amplitude 
and  period  measurements  used  to  compute  station  rn^’s,  more  first  motion 
measurements,  and  more  secondary  phase  data. 

The  question  of  what  improvement  better  or  different  instrumentation 
for  the  stations  of  this  network  may  have  on  the  detection  capability  of 
the  network  was  not  addressed  by  this  report. 
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